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Abstract 

A generalised canonical formulation of gravity is devised for foliations of 
spacetime with codimension n > 1. The new formalism retains n-dimensional 
covariance and is especially suited to 2 + 2 decompositions of spacetime. It is 
also possible to use the generalised formalism to obtain boundary contribu- 
tions to the 3 + 1 Hamiltonian. 
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I. INTRODUCTION 



The canonical methods introduced into the theory of General Relativity by Arnowitt, 
Deser and Misner jl| involve a certain point of view in which spacetime is decomposed into 
spacelike hypersurfaces of constant time. However, there are some situations were a different 
decomposition of spacetime may be prefered and for these situations a new form of canonical 
method is required. 

An important example is provided by the 2 + 2 decompositions of General Relativity [0-[5| 
that arise naturally in the characteristic initial value problem for the Einstein equations. 
Constructing the Hamiltonian along conventional lines requires singling out one of the two 
directions || . It would be desirable to retain some advantages of the Hamiltonian formalism 
without breaking the two-dimensional covariance. 

Another situation of interest arises from considering the thermodynamic properties of 
black holes [[7|. In the usual setup, one envisages a black hole and thermal radiation in 
equilibrium inside a box. The geometry on the surface of the box is fixed when evaluating 
the partition function H8HlQfl. This introduces boundary terms into the Hamiltonian. The 
origin of these boundary terms can be seen in a 2 + 2 decomposition close to the boundary 
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We will describe a generalised canonical method that is appropriate to these situations. 
In the standard canonical formalism, time is introduced by foliating spacetime with spacelike 
hypersurfaces. We will base our canonical formalism on a foliation of codimension n, where 
n < d, the dimension of spacetime. This gives us n coordinates x 1 , not all of them timelike. 

For each variable q l we define n momenta p ! { in terms of the Lagrangian density 
L(q,djq,d a q), where dj are normal derivatives and d a tangential derivatives. We set 

i dL 

and take the definition of the Hamiltonian density to be 

H = p\d I q l -L (2) 

A Legendre transformation can be used to express the field equations obtained from the 
Lagrangian in Hamiltonian form 

V, = f • (4) 

These equations can also be obtained by independent variation of the variables and the 
momenta in the action. 

The case n = d has a history stretching back to De Donder and Weyl |T7|]. This has 
been extended to General Relativity comparatively recently |l8j. The geometrical structure 
of the Weyl-De Donder theory has also been investigated |TP|pO| . We will work within the 
range 1 < n < d, intermediate between the usual canonical approach and the Weyl-De 
Donder theory. 
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The most significant advantage to be gained by introducing a Hamiltonian is the exis- 



tance of canonical transformations. These exist also for the generalised formalism fll9|j20 
and this gives us the freedom to perform changes of variables that simplify the resulting 
field equations. An example of this is provided in section 4, where we consider the 2 + 2 
decomposition of the Einstein equations. 

There are also important limitations to the method. The Hamiltonian density, integrated 
over a hypersurface, does not generate translations along the x 1 directions in the way that 
the usual Hamiltonian does. This complicates the generalisation of Hamilton-Jacobi theory 
and implies that the generalised Hamiltonian density does not have a direct application to 
quantum theory. 

A positive feature of the more general approach is that the momenta have a natural 
variational interpretation. In the usual canonical formalism, the variation of Hamilton's 
principal function with respect to the fields at the final time gives the momenta. In section 
5 we see that, in the new approach, the momenta can be obtained from variation of a 
principal function with respect to fields fixed on a corresponding boundary. This makes the 
new approach a natural way of analysing the partition function for black holes. 



II. LAGRANGIAN DECOMPOSITON 

We suppose that the original manifold M. can be decomposed in n < d coordinate 
surfaces, with coordinates x 1 , and m dimensional subspaces S. Our aim is to decompose 
the Einstein-Hilbert Lagrangian for the metric g, 

L = 7^R\9\ 1, \ (5) 

where k 2 = 8ttG and \g\ = | det g a b\, into a kinetic term, which has derivatives with respect 
to x 1 , and a potential term. 

The decomposed metric g can be written 

g = a ab uj a <g> u b + riuW <g> n J , (6) 

where u) a belongs to the cotangent space of S and n 1 to the orthogonal space. We choose a 
coordinate basis in which n 1 and the dual vectors n/ are given by 

n 7 = dx 1 , nj = 9 Z - Nj a d a , (7) 

where the d a are tangential derivatives. 

When a covariant expression is wanted we replace the coordinate derivative by the Lie 
derivatives Cj along iij. The tangential components of the Lie derivative will be denoted by 
D, 

DiT h - c d ,„ e = ay . . . a c ?af . . . a e e CiT b " c rf /... e / (8) 



The derivative D is the same as the one used in references |2]J and ||. We also define 
tangential vector fields 



[Di,Dj], (9) 
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which depend only on the Nj a . 

Decompositions of the curvature tensor are described in the appendix. These allow us 
to express the Ricci curvature scalar as 

R = 11 R + ± R + kik 1 - k Iab k Iab 

W Ia a J j a -a IJa aj Ia + V a v a . (10) 

where 

k T ab = \rf 3 Djo ah (11) 



a 
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IJ a = \F IJ a + \d a (r l IJ ). (12) 



Also, 

v a = -Infk 1 + 2a 7 /. (13) 

contributes boundary terms. 

The Ricci curvature scalar R can also be expressed in terms of derivatives, 

± R = c IJM c JIM - c\ M c J MJ - 2c 1 {IJ] k J - 2a IJ a a [IJ]a + V b u h (14) 

where 

cijk = -\Di7]jk - \D K r]u + \DjT]ik (15) 

and 

u b = r ] IK {p I r ] j K - D jVlK )n Jb . (16) 

contributes boundary terms. 

Putting all this into the Einstein-Hilbert Lagrangian, we get 

L = T-V (17) 

where 

T = iG^W'D^Dja^ + \G^ K ^ N D Kmj D NVLM + G^ KL ^ b D KVlJ D L a ab 
' 'g\ l ' % a ab F IJ a F IJb (18) 



8k 2 



V = -\G {1J){KL) o ah d amj d mKL - ^- 2 \g\ l ' 2 \\R. (19) 



Gabcd |_|l/2 ( „ab cd „acbd 

= - — -\g\ ' \<7 a —a a 



qIJKLMN = JLUV2 UL JN KM + I[J L]M K N + IJ 
9 k 2 V 



V MN V KL 



The tensors appearing in the kinetic term are 

4/? 
1 

2^ 

G IJKL = ^\g\ 1/2 (v IJ v KL - v IL v JK ) (20) 

symmetrised as indicated by brackets in equation ( fl8|) . 

Similar results to this have been obtained before. The case n = 2 has been analysed by 
Torre 0. Results for general n are quoted by Cho et al. [21], although their Lagrangian still 
retains terms that are second order in derivatives. 
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III. HAMILTONIAN DECOMPOSITION 



The Hamiltonian formalism that we develop here is a generalisation of the usual approach, 
with special attention to retaining the covariance along the tangential and normal directions. 
Momenta are defined in each of the normal directions x 1 conjugate to the quantities a ab , rjjj 
and as 

^ = (21) 
^ = TO 5 ' (22) 

p " = mwf- (23) 

We shall see shortly that these relations can be inverted when the codimension n > 1, 
enabling us to replace derivatives of fields along the normal directions with momenta. The 
case n = 1 leads to the usual Hamiltonian formalism with constraints. 

We define our Hamiltonian density to be 

H = n IJK d Kmj + p Iab d!a ab + P IJ a djNj a - L. (24) 
This can be written in covariant form using the derivative D, 

H = n IJK D KVl j + p Iab Dja ab + \P IJ ' a F LJ a + iV/if^ - L (25) 

where 

H* a = 7T JKI d aVjK - 2aJV cP Ibc + a^F 1 ^^ + P IJ b d a Nj b (26) 

The tangential derivative, defined in the appendix, contains a connection coefficient for the 
normal index and the Fjj a that appears in H I a precisely cancels this coefficient. 

When the derivatives are replaced with the momenta the Hamiltonian density becomes 

H = T + V + N^H^ (27) 

where V is given by equation (|19|) . 



T = \g {ij)k{lm)n ii ijk it lmn + G {/ j) ( ^L)A Ji " L + \G {ab ) { c d) mjP Iab P Jcd 

-l^\g\-^a ab Pu a P IJb (28) 

and p 1 = a ab p ab . The new coefficients are obtained by inverting the matrix 



/ q(IJ)K(LM)N Q{IJ)KL a cd 
y Q(LM)NI a ab Q{ab){cd)^IL 



(29) 



The result of the matrix inversion is that 
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^ 211-1/2/ 2{m — 1) 

Guklmn = 2k \g\ 1 -— — — -tVuVlmVnk + ZvilVjnVkm 

V (d-2)(n-l) 

(VilVjkVmn + VlmVinVjk + VuVklVmn) J (30) 



fn - 1 



G IJKL = 2K 2 \g\' 1 ' 2 [jdZqnJK-niL) , (31) 
Gabcd = -2K 2 \g\~ 1/2 ^2a ac a bd - J^T^) ah ° ' ( 32 ) 

The explicit expressions show that the inversion is possible for codimension 1. When 
n=l, there are m constraints it IJK = P IJ a = 

The Hamiltonian allows us to construct the action in Hamiltonian form 

S = Jd n x (y^dxTju + v lah dia ah + P IJ a djNj a - h) . (33) 

The extrema of this action with respect to independent variations of the momenta 

P IJ a and the fields a a b , ViJi Ni a gives a first order set of equations. The only constraints 

on the momenta are the index symmetries, 

n [IJ\K = p I[ab] = p(IJ) a = q (34) 

If we vary the momenta P IJ a then we find that 

DkVU — G(U)K(LM)N^ LMN + G(u)(KL)P (35) 
Djaab = VijG i a b){cd) p Jcd + G {lm){ni) jc lmn (36) 

F IJ a = 2K 2 \g\- 1,2 P U a (37) 

The covariant derivative appears here because the variation of the H I a term in the Hamil- 
tonian provides the necessary connection components. These equations are equivalent to 
equations ( PH - P^ ) which define the momenta. 

If we vary the action with respect to the fields and substitute for the momenta we 
obtain exactly the same equations as those obtained from the Lagrangian. Therefore the 
extrema of the action give a first order form for the Einstein equations decomposed in n 
independent directions. This is also true in the degenerate case n — 1, provided that we set 

7T IJK EE P IJ n EE 0. 



IV. M + 2 REDUCTIONS 

When the space of normals to the surface is two dimensional then some simplification 
occurs. It even becomes possible to fix some of the gauge freedom and still obtain the full 
set of Einstein equations from the Hamiltonian. 

We replace rjjj by new variables (rj, fju), where 

77 =(1/2) log H (38) 
fju = \v\~ 1/2 Vu (39) 
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This is a canonical transformation if we take conjugate momenta 

vr 7 = VjkK JKI (40) 
~ijk = |^|i/2 (^ijk _ l^W) (41) 

The only term that depends on n IJK in the Hamiltonian has the form 

1 n -IJK-LMN ( AO \ 



where 



Variation of ir IJK leads to 

DkVij = \vr 1/2 GijKLMNir LMN (44) 
If we choose a gauge in which fju is constant initially, then 7t IJK = 0, or equivalently 

tfJK = l^^IJ^K _ ^IK^J _ v JK n Iy (45) 

Variation of fju can be used to show that D k tt ijk = 0, and fju remains constant. Varia- 
tions of the terms involving d a fju and fr IJK vanish and we can drop these terms from the 
Hamiltonian. 

The Hamiltonian becomes 

a tjI 



H = T + V + Nj a H a (46) 



where 



T = 2^\g\- l/2 riu ( — *rV + -pV - a ac a bdP Iab p Jcd + -pV 
\ m m m 

_l K 2 lgl -l/2 a a br]iKT]jL pIJ a pKL b (47) 

V = -^\9\ 1/2 cr ab (daVdbr, + 2WR ab ) (48) 

and we set \g\ = \a\ exp(2rj) and r\u = fju exp(^). 
Variation of the momenta gives the equations 

Djri = ^ 2 \g\- 1/2 vu (~P J + —irA (49) 
\m m / 

D l( r ab = -AK 2 \g\- 1/2 r] I j (p J ab - —a ab p J - —a ab tt j ) (50) 

F IJ a = 2K 2 \g\-^P IJ a (51) 

where Fjj a depends on derivatives of Nj a as defined by equation (|9p. Normal indices are 
raised by r] IJ . 

Variation of the fields gives the equations 
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= ^M 1 / 2 (\\R - "V 2 ^ - 1(11 v v y + ~ w«) 



(52) 



Dj(a ac p Ibc ) = — 



- — 7T J 7r J + — p 7 7T J - /"V^ + — j/j/ 

mm m 



f^M 1/2 ("VJ 1 - c^V 2 ^ + |HV^I!VS - ^("V^) 2 
-^l^| 1/2 (F IJa F IJb + \o a b F IJc F^ 



2k 2 



1 1 a ^ « 



(53) 
(54) 



These equations are far simpler than those derived from the action (33). 

These equations can be put into second order form by eliminating the momenta and 
producing equations for the metric components. Expression for the momenta follow from 



equations (03-51 



V 



lab 



iu|V2 ( k Iab_ a ai k I + l _ a ab D l\ 

Ik 2 \ 2 J 



2k 2 m 



(55) 
(56) 



where the extrinsic curvature ki a b = (^/2)Dja a b. If we differentiate equations (|50"D and (^ 



D^t^D 1 ^ = -2 (Dja acP Ibc )^ + -a b D^ 

\ /TF m 

D^^D 1 ^ = —D^ +(—-2] Djir 1 



m 



m 



where Qtf denotes the traceless part of the tensor. Now we use equations (|53| ) and 
get 

\g\-^Dj (|*| 1/a *V) = ll R a b ~ 11 Vj V b V - ^V a ^W\ + \F IJa F IJb 



\v\ 



■WDj (ivl^D^) =-WR- i(« Vr?) 2 - -F IJa F IJa - (k Iab k Iab - tfkj 



In addition to these equations we still have equation fl54|) , 

Wl^Dj (M^F/.) = (2T7/»V 6 - F/ b )(k Ja b - a b kj - \a b D jr] ) + kjd^ 



(57) 
(58) 

3) to 

(59) 
(60) 

(61) 



Equations (]5"9H6"T|) are a complete set of vacuum Einstein equations and are equivalent to 
the vanishing of the Ricci tensor. They are second order in derivatives of the metric when 
used in conjunction with equations ( |A32| ) for k T ab and (|A31[) for Fjf . 

A similar set of equations, for the 2 + 2 decomposition with null coordinates, has been 
given by Brady et al. H. In their case 



ViJ = e 



-1 
-1 



(62) 
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With the change of notation described in Table 1 of the appendix, and some trivial rear- 
rangements, the two sets of equations agree. One advantage of our more general formalism 
is that we have been able to derive all of the their equations from an action principal. 

Another important example is when the submanifolds have m commuting killing vectors 
and the metric decomposition is orthogonal, Fjj a = 0. In this case, equation ( |6"TD vanishes 
identically and the equations separate, making it possible to solve fl59|) for a a b- The solution 
can then be substituted into equation ( pOD to obtain rj. This special case includes the Kerr 
metric and its generalisations. 



V. BOUNDARY CONTRIBUTIONS TO THE HAMILTON! AN 

We turn now to the application of the canonical formalism to manifolds with boundaries, 
where the metric on the boundary is specified and the boundary terms that were discarded 
earlier will have to be included. This can be readily adapted to the situation in which 
only the intrinsic metric on the boundary is fixed, and we shall see how the action can be 
recovered. 

Problems of this type are often encountered in the path integral approach to quantum 
gravity. The 3 + 1 formalism can be used to decompose the action for manifolds with 
boundary and it has been analysed in this context by York and Hawking and Hunter 



|T5| . We will see how the generalised canonical approach leads to a new understanding of 
the 3 + 1 case and some improvements on the previous results. Because the results have 
applications to quantum gravity, we will make explicit use of Riemannian signature (+ +++) 
in this section. 

The decomposition of the action lead to two total divergences, equations ( jl3|) and (|16|). 



We shall add a boundary term to the original action to cancel these divergences 

S m + n = J_r dR+ J_f dui f 2k I _ IL JK( D _ DjJ]lK j\ ( 63 ) 

2k z Jm 2k z JdM v ' 

This assumes that the boundary is given by a function of the x 1 only. We set duji = nfduja, 
where duj is the surface form and dfi is the volume measure. 

After canonical decomposition the Lagrangian is first order in field derivatives. Variation 
of the action gives the field equations, plus a boundary term 

(5S m+n )dM = I |^r 1/2 ^/ U IJK 5 VJK + p Iab 5a ab + P IJ a 5N/) . (64) 

The field equations can be deduced from the variational principal if the metric is fixed on 
the boundary. We can also take the value of the action for solutions to the field equations, 
calling this the principal function, then equation (|)4]) relates the momenta to the variation 
of the principal function. 

The action S m+n is suitable for fixing all of the metric components on the boundary but 
it is not yet suitable for variations which fix only the intrinsic components of the metric. 
We will obtain the necessary modifications to the action for the boundary shown in figure || 
and compare this with known results. 

We assume that the space can be foliated by two sets of coordinate surfaces which include 
the boundary, allowing a 2 + 2 decomposition. Figure |2| shows part of the boundary, where 
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n 7 are normal forms and the dual vectors n/ are tangential to the boundary surfaces. We 
let 9 be the angle between the tangential vectors, then 9 can be related to the metric rjjj by 

V12 = (7711*722) 1/2 cos 9 (65) 
l^| 1/2 = (?7ii^22) 1/2 sm6> (66) 

We also make use of the traced extrinsic curvature of the boundary K 1 = Vn 7 . Comparison 
with the definition (|A3|) gives 

K 1 = k ! + c JI J -(r l II )- l c in (67) 

with no sum over I. 

With a 2 + 2 decomposition it is possible to take r/jj to be conformally constant, and by 
equation fl45]) 



7T 



UK _ ^IJ^K _ l v JK n I _ yK^J (6g) 



This simplifies equation (j64|) , 

t, ijk 5t 1ij = (r] IJ 7T K - rfK^briu (69) 
The antisymmetry in the indices / and K allows us to introduce the antisymmetric tensor 

n IJK 5 mj = e IK e JL ir L 6 VlJ (70) 
The momenta can also be replaced by the extrinsic curvatures 

^ JK Svu = ^ 1 \g\ 1/3 ^ L k L Smj (7i) 

Since this vanishes when I = K, the variation of the action does not depend upon the purely 
normal metric components but it does depend upon the mixed components #7712- We can 
eliminate these by using a new action S cov , obtained by adding an extra surface term. 

Only one extra surface term, depending on 9 and kj, is needed. The surface term becomes 
(compare equation (|63|)) 

/ dcoj (2k 1 - V IL V JK (DjVlk ~ DjVik) - 0e IJ kj) (72) 

ZK Z JdM v J 

Variation of the final term cancels the br\\i- 

We can rewrite the boundary term in terms of the extrinsic curvature K 1 of the boundary 
using the definition flB"7p. We also use the expression ( |A35| ) for cjjk, to get a familiar 
expression 

SC ° V = I d ^ R + \ f dujjK 1 + \ [ dti(0 o - 9) (73) 

2fiT JM K JdM K JT, 

where S denotes the corners of the boundary, d[i the volume measure and 9 is a constant. 
For an additive action we take 9q = ir/2. 
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The corner terms were first discovered in the context of Regge calculus by Hartle and 
Sorkin p2| . They were rediscovered later with Lorentzian signature [13,14]. In the Lorenzian 
case we replace 9 — 9 by i/3, where j3 is a Lorentz boost. 

For the remainder of this section we will decompose S cov into 3 + 1 form and construct 
the corresponding Hamiltonian. We will take x 1 to be the time coordinate, running over a 
fixed interval, and x 2 a radial coordinate, running from zero to a fixed value x 2 s on the edge 
S. The coordinates on the surfaces of constant time will be denoted by x l , and 

x 2 (x\t) = x 2 s (74) 

on S, as shown in figure (|3]). The metric becomes 

g = iV 2 dt <g> dt + ^-(dx* + N*dt) <g> (dx. j + N j dt). (75) 

The 3 + 1 and the 2 + 2 form of the metric are related by 77 11 = N~ 2 and iV la = h ai N\ 

The first step is to construct the Lagrangian. We separate the volume terms from the 
boundary terms and write S cov = Sm + SftM- The 3 + 1 decomposition of the curvature 
introduces divergences V • v, with v given by of equation (|13D in codimension 1. This 
introduces boundary terms of the form n / • v, which become 

n 1 ■ v = -2K 1 (76) 
n 2 -v = 2(r / 11 )- 1 ((^ 2 )- 1 c 11 2 -irV 2 ) (77) 



Boundary terms on the surfaces of constant time cancel with boundary terms already present 
in the action S cov (equation 173]). Combining terms on the boundary S with S cov gives 

SdM = h J/^T'iM^c^-KW' + V 11 ^ 2 ) (78) 
There are also corner terms. Next we replace the extrinsic curvatures using the identity 

k 2 = V21K 1 - V n \ V \ 1 / 2 D 1 9 + W^lc 11 * (79) 
which can be derived from from equation (p7|) and ( |A35| ) . We are left with 

S 9M = ^J s dfxdt ([r 1 ll r l22 )- l ' 2 k 2 -(9- 9 )h) (80) 

When #0 = 7r/2, the corner terms cancel and this expression is the complete expression for 
boundary contributions to the Lagrangian. 

The surface term ( |S0|) has two components. Appart from at the corners, the value of k\ 
is fixed by the boundary conditions. However, k\ includes time derivatives of the metric, 

h = \a ab (d t a ah - 2"V a iVi & ) (81) 

When we take the functional derivative of the the Lagrangian with respect to d t u a b we get 
surface momenta, 

P ab = ^W\ 1/2 (9 -9)a ab (82) 
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This suggests that the angle 9 should be regarded as a momentum variable. Since the value 
of 9 on the corners depends upon the normal metric components it is not fixed by the 
boundary conditions. 

In canonical form, the surface terms (|80|) become 

S 9M = I d 2 xdt (v ah d t a ab + Nk- 2 ^)- 1 / 2 ^ 2 ^ + N l hJV b p ab ) (83) 

The usual ADM Hamiltonian can be obtained by specialising the earlier discussion ton = 1, 
setting tt ijk = 0. This introduces an additional surface term from a divergence 

- 2< 3 >V i (iV i p y ) (84) 

Combining the surface terms together gives 

/ d 2 xdt(p ab d t a ab - NH- WHi) (85) 

The surface contributions to the Hamiltonian are, 

n = -K- 2 { n22 y l i 2 D 2 \a\ 1/2 (86) 

Hi = -hJV b p ab - (r/ 22 )- 1/2 ^n 2 ^ fc (87) 



Variation of the full action, subject to the condition (|8"2"|), gives the field equations with no 
surface terms. 

The term Tli looks very different from the results of Hawking and Hunter |15] because 
we have separated the momentum term in (|83|). On the other hand, the value of the action 
is identical in the two cases. 

We can further verify the existence of surface momenta by direct variation of the action 
S cov . This gives the field equations plus surface terms 

SS cov = I (PxpVSha + I d 2 xp ab 5a ab . (88) 

JdM JT, 

The value of S cov for solutions to the field equations gives Hamilton's principal function. This 
will be a function of the metric on DM.. Equation fl64|) shows that the principal function has 
a variation with support on the corners. Hamilton- Jacobi theory tells us that this variation 
has to be a momentum. It also follows that one family of constants of integration of the 
Hamilton-Jacobi equation will be the initial values of the surface momenta, or equivalently, 
the initial values of 9. 



VI. CONCLUSION 

We have described a generalisation of the Hamiltonian formalism for gravity that is 
appropriate to foliations of codimension n > 1. The new formalism puts the vacuum Einstein 
equations into a first-order form. It also has a variational formulation which allows canonical 
transformations to be used to simplify the classical field equations. 

The results presented here can be extended to include matter fields. We hope that they 
will be valuable for finding new solutions to Einstein-matter systems in various spacetime 
dimensions. 
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APPENDIX: GENERAL FORMALISM 



Suppose we have a metric g on a d-dimensional spacetime M. and a sub-surface E of 
dimension m < d. Using a set of normal forms n 7 and tangential forms u> a we can put the 
metric in the form 

g = a ab LJ a <g> u) b + riuii 1 <g> n J , (Al) 

The Levi-Cevita connection V on M. induces a connection on the tangent space to £ and 
we can obtain Gauss-Codazzi relations for the curvature. There does not exist a standard 
notation for these reductions and so we set out ours here. 
The inverse metric 

g = a ab e a <g> e b + ?/ J n 7 <g> nj, (A2) 

where the vectors e a and nj are dual to the forms u> a and n 7 . We call tensors with compo- 
nents along e a or u> a tangential and those with components along n/ and n 7 normal. 
We define the extrinsic curvature and related quantities, 

k I ab = a a a 'a b b '(V a/ n I ) b ,, (A3) 

a IJ a = a a a 'n Ib '(V a ,n I ) b ,, (A4) 

b IJ a = a a a 'n Ih \V al W) bl) (A5) 

c IJK = n Ka 'n Ib '(V a ,n I ) b , (A6) 

These transform like connections under changes in the normal forms, but they can be re- 
garded as rank 3 tensors if the basis is changed whilst keeping the normal forms fixed. 

The hypersurface orthogonality of the n 7 implies that the extrinsic curvature is symmetric 
in the tangential indices. The tensors b IJ a and c IJK are identically zero in the most familiar 
situation where the surface has codimension 1 and the normals have unit length. In the 
general case they can be interpreted as connection coefficients in the normal directions. 
This can be seen in the Gauss- Weingarten equations that follow from the definitions given 
above, 

Vn 7 = k 7 + a/ b n J <g> u) b + b/ a uj a <g> n J + c] K n J <g> n K (A7) 
Vn 7 = k/ + aj/ V ® e b - b/ a uj a <g> nj - c/ jn J <g> n K (A8) 

We can now decompose the connection in directions tangential and normal to the surface, 

^ a T J '" K L...M b "' C d...e = ^ a'T J "' K L...M h "' C d...e (A9) 

V/V 7 -'^/ ,/••',,., = nr'v./r'^, ,/ (Aio) 

The derivatives have been projected so that tangential indices remain tangential and normal 
indices remain normal. These decomposed derivatives are both metric connections 

"V a (7 6c = ^itjjk = 0. (All) 
However, the perpendicular connection has torsion a' 77 ^. 
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With these derivatives, we can now complete the Gauss-Weingarten equations and de- 
duce that 

Vn 7 = "Vn; + ^Vnj + k Ia b u a <g> e b + ajI b n J <g> e b (A12) 
Ve 6 = 11 Ve 6 + L Ve b - k 1 ab u a <g> n 7 - a/ fe n 7 <g> n, 7 (A13) 

It is often convenient to replace -'-V and use instead the tangential projection of the Lie 
derivative along n/, which we denote by Dj. The vanishing torsion of the Levi-Cevita 
connection V implies that the two are related by 

^Vefe = De b + k Ib c n T <g> e c (A14) 
= Du) b - k^n 1 <g> u) c (A15) 

The curvature tensors of the connections will be denoted by "R and -"-R respectively. 
The perpendicular derivative of a normal vector X is given by 

^V/X 7 = n/ (X J ) - c/ K X K (A16) 

The commutator of two of these derivatives gives an expression for -'-R in terms of the 
Gauss- Wiengarten coefficients, 

L R IJ KL = 2 ± V [L c x] / J + 2 C[i f J c L \ M + 2b JI a a [KL] a (A17) 

The final term, that depends on the torsion, comes from the commutator [n^-, nj. Similarly, 
we also have an expression for the commutator of two tangential derivatives "V acting on 
normal vectors, 

ll R IJ ab = 2^ [a b IJ b] - b IK a b J Kb + b IK b b J Ka (A18) 

Taking commutators of V on the Gauss-Weingarten relations we find Gauss-Codazzi 
relations for the tangential a,b . . . and normal components I, J . . . of the curvature tensor 
R of the connection V, 

R\ cd = ll R\ cd - Vu (k Ia ck J bd - k Ia d k\ c ) (A19) 

Rjbc = "V^'ac - ll V c k T ab (A20) 

R T a \ = -^ J k ! ab - k 1 ^^ + WV b a JI a - a K \a JK b . (A21) 

R IJ ab = ll R U ab ~ k Ic a k J bc + k Ic b k J ac (A22) 

R a UK = ± V J a X / a - "V^aj/ - 2kj a b a [JK] b (A23) 

R 1 jkl = ±rI jkl - a Kb a LJ b + a L \a K b . (A24) 

In our notation, the covariant derivatives include connection coefficients for the normal 
indices. 

By taking traces of the above relations, we have a decomposition of the Ricci scalar R: 

R = \\R + ± R + kik 1 - k Iab k Iab 

W Ia a J j a - a IJa a JIa + V a v a . (A25) 



14 



where 

v a = -In^k 1 + 2a 7 /. (A26) 

A special case of interest is when the manifold M. has a foliation of codimension n. In 
this case we can choose a coordinate basis, 

n 1 = dx 1 , u a = dx a + Nfdx 1 . (A27) 

The dual basis 

n T = d T - Nj a d a , e a = d a , (A28) 

where dj = d/dx 1 and d a = d/dx a . 

In this basis the Gauss-Weingarten coefficients are related by a IJ a = b IJ a and c IJK = 
c KIJ . Because of these relationships the derivative Dj satisfies 

D/n/ = D^ 1 = 0. (A29) 

Therefore there are no connection terms for the normal indices. Furthermore, the tensor 
Fjj, defined by 

Fu = [Dj, Dj] (A30) 

can also be written 

Fu a = djNj a - djNj" - Nj b d b Nj a + Nj b d b Nj a (A31) 

It depends only on the A^ a . 

Now we can express the Gauss-Weingarten coefficients in terms of the derivative Dj. 
First of all, using ^Vaab = and equation (|A15| ), we find that 

k T ah = lv IJ Dja ab (A32) 

Writing out the Lie derivative in the coordinate basis gives, 

Dja ab = m (a ab ) - (T ac d b Nf - a bc d a Nf (A33) 

In the coordinate basis we get 

a IJ a = b IJ a = lF IJ a + ±d a (r ] IJ ). (A34) 

Finally, by regarding the c IJK as connection coefficients we may deduce that 

cijk = -\F>iT]jk ~ \Dkthj + \Dji] IK . (A35) 

The table gives comparisons between the notation used here and the notation used in two 
other references. 
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TABLES 



GM BDIM HH 



-e 



/ —1 sinhr? 



\ — e x ,/ Vsinhr/ ! 

Cab 9ab Cab 

n/ I a Ti sechr/ 

n 1 e~ x £ A n M 

ab K^ab ^ab 

^V a V a (see caption) 

Di D A 
au a L AB a 

CIJK N B CA 

TABLE I. This table compares the notation of the present paper (GM) with reference (BDIM) 
and (HH). The index a = 1 . . . m, where BDIM and HH both have m = 2. The derivative " V a 
includes connections on normal indices but V a does not. 
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FIGURES 



x 2 const 




x 1 const 



£ 

FIG. 1. The manifold M with a boundary made up from two surfaces of constant x 1 and two 
surfaces of constant x 2 intersecting in four corners of dimension two. 




FIG. 2. Surfaces of constant x 1 , I = 1,2, with the basis of normal vectors n 1 and the dual 
basis n/. The angle between the dual basis vectors is 6 




FIG. 3. The manifold M. with boundary S and two surfaces of constant time. 
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